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Abstract: Water-use efficiency (WUE) is a key plant functional trait that plays a central role in the global 
cycles of water and catbon. Although increasing precipitation may cause vegetation changes, few studies 
have explored the linkage between alteration in vegetation and WUE. Here, we analyzed the responses of 
leaf WUE, ecosystem carbon and water exchanges, ecosystem WUE, and plant community composition 
changes under normal conditions and also under extra 15% or 30% increases in annual precipitation in a 
temperate desert ecosystem of Xinjiang, China. We found that leaf WUE and ecosystem WUE showed 
inconsistent responses to increasing precipitation. Leaf WUE consistently decreased as precipitation 
increased. By contrast, the responses of the ecosystem WUE to increasing precipitation are different in 
different precipitation regimes: increasing by 33.9% in the wet year (i.e., the normal precipitation years) 
and decreasing by 4.1% in the dry year when the precipitation was about 30% less than that in the wet year. 
We systematically assessed the herbaceous community dynamics, community composition, and vegetation 
coverage to explain the responses of ecosystem WUE, and found that the between-year discrepancy in 
ecosystem WUE was consistent with the extent to which plant biomass was stimulated by the increase in 
precipitation. Although there was no change in the relative significance of ephemerals in the plant 
community, its greater overall plant biomass drove an increased ecosystem WUE under the conditions of 
increasing precipitation in 2011. Howevet, the slight increase in plant biomass exerted no significant effect 
on ecosystem WUE in 2012. Our findings suggest that an alteration in the dominant species in this plant 
community can induce a shift in the carbon- and water-based economics of desert ecosystems. 


Keywords: desert ecosystem; ecosystem water-use efficiency; gross carbon exchange; increasing precipitation; leaf 
water-use efficiency; net carbon exchange; Gurbantunggut Desert 


1 Introduction 


Precipitation is the determining factor for plant growth, community composition, and ecosystem 
functioning in deserts (Webb et al., 1983; Austin et al., 2004; Lal, 2004). Changes in the 
precipitation amount, frequency, interval, and seasonality can induce substantial alterations in the 
carbon cycling and water use strategies of plants (Robertson et al., 2010). However, it remains 
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unclear how variation in precipitation affects the carbon exchange and water-use efficiency (WUE) 
at the ecosystem level (Huxman et al., 2004; Potts et al., 2006; Thomey et al., 2011). 

WUE is a comprehensive index of carbon and water fluxes and is recognized as an important 
functional trait of plant individuals and communities (Le Houerou, 1984; Linderson et al., 2012). 
WUE is closely related to plant photosynthesis and evapotranspiration, and thus represents a plant 
adaptive strategy at both the leaf and ecosystem scales (Eamus et al., 2013). Changes in WUE at 
different scales may reflect the responses of plants to different ecological processes. Leaf WUE is 
mainly controlled by plant stomata (Toft et al., 1989; Golluscio and Oesterheld, 2007), and leaf 
WUE responses to increasing precipitation depend on the relative changes in photosynthesis and 
transpiration (Cohen, 1970; Golluscio and Oesterheld, 2007). Ecosystem WUE integrates the 
comprehensive responses of vegetation composition and ecosystem evaporation for a given 
community (Law et al., 2002). Although ecosystem productivity is positively correlated with 
precipitation across a sufficiently large spatial scale, ecosystem WUE may either increase or 
decrease with increasing precipitation (Varnamkhasti et al., 1995; Veron et al., 2002; Scanlon and 
Albertson, 2004; Hastings et al., 2005), and the consensus on how ecosystem WUE responds to 
increasing precipitation is still lacking (Le Houerou, 1984; Bell et al., 2012; Liu et al., 2012). The 
lack of consensus arises probably because the photosynthetic assimilation of dominant plant 
species often determines the magnitude of net primary production (NPP) responses to increasing 
precipitation (Golluscio and Oesterheld, 2007). Moreover, vegetation coverage further influences 
the proportion of plant transpiration that contributes to ecosystem evaporation (Huxman et al., 
2004). Since water and carbon fluxes are coupled at the ecosystem level and both depend on plant 
community composition, studies on the variation of WUE at different scales (leaf vs. community) 
can reveal the internal responses of ecological process to background climate changes (Linderson 
et al., 2012). 

The Gurbantunggut Desert in northwestern China represents an important part of the temperate 
desert biome in Central Asia, one that is predicted to experience increasing precipitation from 
2001 to 2030 (Liu et al., 2010). The objectives of this study were threefold: (1) to test the 
responses of WUE, at both leaf and ecosystem levels, to increasing precipitation and to analyze 
the relationships between WUE at the two scales; (2) to explore how community composition 
changes under increasing precipitation might influence ecosystem WUE; and (3) to test whether 
the WUE of the herbaceous community under shrub canopies and in interplant habitats exhibits a 
consistent response to increasing precipitation. Specifically, we tested two proposed hypotheses: 
first, ephemerals and summer-annuals are the two major plant functional groups in the 
Gurbantunggut Desert, which undergo a community shift in their dominance from ephemerals in 
April-early June to summer-annuals in early June-August. Since the relative density of these 
plants and their phenology also change with increasing precipitation, we hypothesized that WUE 
at both the leaf and ecosystem levels would increase in response to increasing precipitation for 
sustaining ephemeral growth. The WUE at community level, however, would decrease the limited 
stimulation of summer-annual growth due to low productivity and high soil respiration and 
evaporation. Second, our study site exhibits a typical mosaic landscape, with “fertilized islands” 
occurring under shrub canopies. The canopies are featured by cyanobacterial-lichen crusts that in 
turn cover the soil in the interplant matrix (Su et al., 2007, 2013; Li, 2012). Since plant biomass in 
the interplant space was higher than that found beneath the shrubs in a normal precipitation year, 
but similar in a dry year (Huang et al., 2015), we hypothesized that ecosystem WUE for the 
interplant area would also exceed that of the beneath-shrubs area during a normal precipitation 
year, while both (interplant area and beneath-shrubs area) are similar in a dry year. 


2 Materials and methods 


2.1 Study area 


This field experiment was conducted at the margin of the Fukang National Field Scientific 
Observation and Research Station for Desert Ecosystems, Chinese Academy of Sciences, on the 
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southern edge of the Gurbantunggut Desert in the hinterland of the Eurasian continent (44?17'N, 
87°56'E; 475 m a.s.l). This region has a continental arid temperate climate, with a hot, dry 
summer and a cold winter. The annual mean temperature is 6.6°C, while the mean annual 
precipitation is 160 mm, of which 70%-—80% occurs during the plant growing season from April 
to September. The multi-year winter average of snow thickness is 20 cm. The soil here is 
classified as chestnut soil (Kastanozem), and is well drained with a low-moisture holding capacity. 
There is ca. 30% vegetation cover, with the dominant plant species being Haloxylon 
ammodendron, H. persicum, and Tamarix ramosissima. The herbaceous layer is primarily 
composed of Alyssum linifolium, Astragalus arpilobus, Ceratocarpus arenarius, Descurainia 
sophia, Erodium oxyrrhynchum, Euphorbia turczaninowii, Hyalea pulchella, Lactuca undulate, 
Orostachys spinosus, Schismus arabicus, Trigonella arcuate, and Agriophyllum squarrosum. 
Similar to other desert ecosystems, the precipitation regime dramatically influences primary 
productivity in this desert, being strongly associated with inter-annual and intra-annual 
precipitation patterns (Schimel et al., 1990; Ojima et al., 1998). 


2.2 Experimental design and soil sampling 


We used a completely randomized block design with three treatments: natural precipitation (C), 
an extra 1596 (P15), and an extra 3096 (P30) of precipitation added, with each treatment replicated 
six times. A total of 18 plots were arranged in a 4x6 matrix; each plot was 10 mx10 m, with a 
10-m-wide buffer zone between adjacent plots. For the P15 and P30 treatments, precipitation was 
increased by 15% and 30% because these are predicted for northern China over the next 30 years 
from 2001 onward (Liu et al., 2010). Precipitation was collected in rainfall collection pans. These 
pans were constructed from galvanized iron sheets, each with a "catch" area of 1.5 mx1.0 m. A 
total of 9 and 18 pans were installed in each 100-m? plot of the P15 and P30 treatments, 
amounting to ca. 15% and ca. 30% of the plot area, respectively. The pan was erected at an angle 
of 2 degrees, and the intercepted rainfall was collected in a bucket buried in the soil. Immediately 
after a precipitation event, the collected rain was evenly sprayed onto the plots. Further, given the 
importance of snow for the herbaceous community dynamics at our study site (Fan et al., 2012), 
snow that had fallen into the pan was also evenly added to the corresponding plot in early spring. 
The experimental design and instrument arrangements were applied in August 2010. In addition, 
in late March 2011, two square aluminum frames (each 0.5 mx0.5 m) were put into the soil at a 
3-cm depth in each plot: one in the interplant matrix, and the other beneath the shrubs of H. 
ammodendron. 


2.3 Meteorological data, soil moisture, and temperature measurements 


Precipitation and air temperature were monitored hourly by an automatic meteorological station 
(Campbell Science Equipment, Logan, USA). Soil moisture at 0-10 cm depth was measured 
using time domain reflectometry (Diviner-2000, Sentek Pty Ltd., Balmain, Australia), while the 
soil temperature at 0—5 cm depth was measured using a thermocouple device (Microlab Lite, 
Fourier Technologies Ltd., Rosh HaAyin, Israel) Two probes, one for the time domain 
reflectometry and the other for the thermocouple, were installed in each plot, with one in the 
interplant space and the other beneath the shrubs of H. ammodendron. Soil moisture and 
temperature measurements were conducted on the same day as CO2 and H20 flux measurements. 


2.4 Carbon and water exchange measurement, and ecosystem WUE calculation 


Ecosystem CO» and H20 fluxes were measured with a transparent chamber (0.5 mx0.5 mx0.5 m) 
attached to an infrared gas analyzer (IRGA; LI-840, LiCor, Lincoln, USA). This method has been 
widely used to measure ecosystem CO» and H20 exchanges in previous field studies (e.g., Risch 
and Frank, 2007; Xia et al., 2009; Niu et al., 2011). The chamber walls were constructed by clear 
polyethylene sheeting, which allows 90% of photosynthetically active radiation to pass through it. 
Two electric cooling fans were installed to circulate the air in the chamber while the 
measurements were taken. CO» and H2O fluxes were measured at 09:30-11:30, at 2-week 
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intervals, during the plant growing seasons of 2011 and 2012 (1.e., from the start of April to the 
end of September). The chamber was sealed to the frame surface, and the two small electric fans 
ran continuously to promote air mixing within the chamber during and between measurements. 
Twelve consecutive recordings of CO» and H20 were taken, at an interval of 10 s during a 120-s 
period, after steady-state conditions were attained within the chamber. The observed increases in 
air temperatures were usually «0.5?C within the chamber during the period of measurement. The 
CO» and H20 concentrations were allowed to build up or decrease over time; the CO» exchange 
and H20 efflux were determined from the time-course of the concentration to calculate net 
ecosystem exchange (NEE) and evapotranspiration (ET). Following the measurement of NEE, the 
chamber was vented and replaced on each frame and covered with an opaque cloth. Then, the CO2 
exchange measurements were repeated. Because the second set of measurements eliminated light 
transmission (and hence photosynthesis was inhibited), these values obtained represented 
ecosystem respiration (ER), which is composed of plant dark respiration and soil respiration. The 
difference between the NEE and ER values was considered to represent the gross ecosystem 
exchange (GEE) for the vegetation within the chamber. NEE, ER, and ET were calculated 
following the formula below, described in the LI-840 manual: 

__10VP,I-W,)_ ôC (1) 

RST, +273.15) ôt’? 

where F. is the CO» (umol/(m?-s)) or H2O (mmol/(m?.s)) exchange rate; V, the volume of the 
assimilation chamber (75,000 cm?); Po, the initial pressure (kPa); Wo, the initial water vapor mole 
fraction (umol/mol); R, the universal gas constant (8.312 kJ/(mol-K)); S, the soil surface area 


(2500 cm?); To, the initial air temperature (°C); and 2C, the slope of least squares linear 
ôt 


regression of CO2 or H20 concentration as a function of time. 

By convention, the negative CO2 exchange values refer to CO2 release, whereas the positive 
values represent CO2 uptake at the sampling site. Ecosystem WUE was calculated as GEE/ET 
(i.e., WUEcre) and NEE/ET (i.e., WUEnre) in our study. 


2.5 Photosynthesis measurements 


Leaf photosynthesis and transpiration rate measurements were only conducted for three dominant 
species, including the two spring ephemeral species (E. oxyrrhynchum and A. linifolium) and one 
spring-summer annual (C. arenarius). Growth peaks occurred in the middle of May for the spring 
ephemeral and in mid-August for the spring-summer annual, so our measurements were taken 
between 10:30 and 11:30 on 15 May, 2011 and 11 May, 2012 for spring ephemerals, and likewise 
on 20 August, 2011 and 2012 for the spring-summer annuals. Leaf photosynthesis (Pn) and 
transpiration rate (TL) were measured with a 6-cm? clamp-on leaf cuvette attached to an open 
gas-exchange system (Li-6400; Li-Cor Inc., Lincoln, USA). Because the leaves of the three herbs 
were very small, typically <1 cm? in area for a single leaf, usually 3—5 leaves were placed into the 
leaf cuvette. In arranging the leaves, they were not overlapped, so as to correctly measure P, and 
Ti. During these measurements, the leaves were illuminated at 900 jumol/(m?.s), using a LED 
light system. Immediately after taking the measurements, the leaves were cut with a pair of 
scissors and placed inside an envelope; samples were then stored at 4°C in a refrigerator for 
subsequent leaf area measurements. We calculated the values of Pn and Ty based on the estimated 
leaf area by taking photos. Leaf WUE was calculated as WUE;zP,/Tz (Jones, 2014). 


2.6 Vegetation biomass and WUE calculation 


Plant biomass was estimated using allometric equations that were obtained by the destructive 
sampling of 750 individuals per species in 2010. The allometric equations for E. oxyrrhynchum, A. 
linifolium, S. arabicus, S. foliosa, and C. arenarius were | TAGB-0.0733e9 1327H. 
TAGB-0.0082e9. 757, TAGB-0.0177e918!7, TAGBZ0.0359e* 7?" and TAGB-0.115e02657H, 
respectively, where TAGB and H express the total plant aboveground biomass (g) and plant height 
(cm), respectively. All five R? values for the allometric equations were 0.87. In 2011, five 1 mx1 m 
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subplots were set up in each plot: two situated beneath shrubs, and three located in the interplant 
matrix space. Herbaceous plant height and individual counts of the five species were measured at 
growth peak in 2011 and 2012, respectively. According to one recent study, the biomass of the 
same five species amounts to more than 85% of the total community biomass at the study site 
(Huang and Li, 2015). 


2.7 Statistical analysis 


A repeated-measure ANOVA was used to test the effects of year (Y), water addition (W), 
microsite (S), and their possible interactions on six response variables: soil moisture, soil 
temperature, GEE, NEE, ET, and ecosystem WUE. If there was a significant inter-annual 
variability (i.e., a year effect, P«0.05), then separate two-way ANOVAs were used to further 
examine the main and interactive effects of water addition and site on the variables in 2011 and 
2012. For Pp, Ti, and leaf WUE, a repeated-measure ANOVA was used to test the effects of Y, W, 
S, and their possible interactions. For both of the two statistics, we initially included block as a 
random factor, it was then removed because it was insignificant in all cases. Pearson correlation 
analysis was used to test the linear associations of WUEcrr, WUEwzr, GEE, NEE, ET, and plant 
biomass with annual precipitation, those of GEE and NEE with ET, and those of Pn, Tr, and 
WUE, with the soil moisture of the three plant species. All the statistical analyses reported here 
were conducted with SPSS software (SPSS for Windows, v13.0, Chicago, USA). 


3 Results 


3. Microclimate 


Total annual precipitation differed between 2011 and 2012, being 167.4 and 102.0 mm, 
respectively. Precipitation in 2011 was rather close to the long-term average annual precipitation, 
while precipitation in 2012 is lower than the average annual precipitation. Thus, we can consider 
the year of 2012 as a dry year. Meanwhile, precipitation in the two years showed strong seasonal 
variations, with 74.476 and 48.896 of the annual precipitation falling during the plant growing 
season in 2011 and 2012, respectively. The average soil volumetric water content (SVWC) in the 
control treatment was 3.44% in 2011 and 3.63% in 2012 (Figs. 1a-d). The P15 treatment did not 
significantly improve SVWC, but the P30 treatment markedly increased SVWC by an average of 
1.9% in 2011 and 2.9% in 2012. Soil temperature (Ts) differed significantly between the two 
years (Figs. le—h; Table 1; P«0.05), but increasing precipitation exerted no discernible influences 
on T; (Figs. le-h; Table 1). 
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Fig. 1 Seasonal patterns of soil volumetric water content (SVWC) at 0-10 cm depth and soil temperature (Ts) at 
0—5 cm depth in the control (C), extra 1596 precipitation addition (P15), and extra 3096 precipitation addition (P30) 
in the interplant matrix and beneath the shrubs of Haloxylon ammodendron (“Beneath H. a") in 2011 and 2012. 
Bars represent standard errors. 
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Table 1 ANOVA results for the effects of year (Y), precipitation (P), microsite (S), and their interactions on soil 
volumetric water content (SVWC), soil temperature (Ts), water-use efficiency (WUEcee and WUEnep), 
herbaceous aboveground biomass (Biomass), gross ecosystem carbon exchange (GEE), net ecosystem carbon 
exchange (NEE), and ecosystem evapotranspiration (ET) 


Effect SVWC Ts WUEcEE WUEnceE Biomass GEE NEE ET 
Y 6.4" 4.5" 37.0 25.67 70.4** 61.3** 75.37 32.5" 
YxP 1.9 23 0.9 5.1° 6.4" 2.5 7:5'* 0.7 
YxS 1.1 2.1 0.0 0.6 2.1 0.0 0.6 0.1 
YxPxS 0.5 2.3 0.0 0.9 0.4 0.9 3.0 0.2 

P 20.4"* 2.4 Oo" 3.4* 20.5"* 10.5** 3.9* 8.7" 
S 0.5 1.7 0.6 3.0 1.2 4.2* 2.3 2.1 
PxS 0.7 0.4 0.9 1.8 2.4 4.8* 4.1* 3.9* 
Note: "^5, **, “denote significance at P<0.001, P<0.01, and P«0.05 levels, respectively. 


3.2 Leaf-level WUE in response to increasing precipitation 


Increasing precipitation significantly increased both P; and TL, notably so for E. oxyrrhynchum (E. 
o) and A. linifolium (A. I) in 2011 and C. arenarius (C. a) in 2012 (Figs. 2a-d). Increasing 
precipitation slightly decreased leaf WUE of all three species, but this was not significant (Figs. 
2e and f). 


P: P=0.02 é P: P=0.4 
S: P«0.001 S: P«0.001 
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Fig. 2. Maximum leaf photosynthesis (Pn; a, b), leaf transpiration (Tr; c, d) and leaf water-use efficiency (WUEL; 
e, f) of two spring ephemerals, E. oxyrrhynchum (E. o) and A. linifolium (A. 1), and one spring-summer annual, C. 
arenarius (C. a), in response to increased precipitation in 2011 and 2012. Measurements were taken on 20 May 
for the spring ephemerals and on 15 August for the spring-summer annual. Bars represent standard errors. 
Different lowercase letters indicate a significant difference at P«0.05 level among precipitation treatments for a 
specific species. C, control; P15, extra 15% precipitation addition; P30, extra 30% precipitation addition. P, 
precipitation; S, species. 


3.3 Ecosystem-level WUE in response to increasing precipitation 


WUEazz and WUEwzz showed strong seasonal variations, ranging from 0 to 1.97 umol CO2/mmol H20 
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for WUEazrz and from —2.37 to 0.31 umol CO2/mmol H20 for WUEneg, respectively (Fig. 3). On 
average, P15 promoted WUEsczz by 13.1%, and P30 promoted WUEazz by 34.1% (Fig. 4). 
Increasing precipitation effects on WUEwzE was dependent on experimental year (Table 1). 
However, WUEnez only significantly increased with increasing precipitation in 2011, on average 
by 18.9% and 33.6% respectively under the P15 and P30 treatments, while no significant impacts 
were found in 2012 (Table 1). WUEgze and WUEmwzz were similar between the microsites of 
interplant matrix and beneath H. ammodendron shrubs (Table 1). 
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Fig. 3 Ecosystem water-use efficiency (WUEcer (GEE/ET, a-d) and WUEnre (NEE/ET, e-h)) in response to 
increasing precipitation in the interplant matrix and beneath the shrubs of Haloxylon ammodendron (“Beneath H. 
a”) in 2011 and 2012. Bars represent standard errors. C, control; P15, extra 15% precipitation addition; P30, extra 
30% precipitation addition. 
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Fig. 4 Effects of increasing precipitation on the ecosystem water-use efficiency (WUEcex (a) and WUENwrz (b)) in 
the interplant matrix and beneath the shrubs of Haloxylon ammodendron (“Beneath H. a") in 2011 and 2012. Bars 
represent standard errors. Different lowercase letters indicate significant differences at P«0.05 level among the 


precipitation treatments within a species. C, control; P15, extra 15% precipitation addition; P30, extra 30% 
precipitation addition. 


3.4 Responses of plant community composition to increasing precipitation and its 
correlation with WUE 


Increasing precipitation slightly increased the species richness of the plant community, but not in 
a Statistically significant way (Table 2). Nevertheless, the densities of the dominant species (E. 
oxyrrhynchum, Alyssum linifolium, S. arabicus, Ceratocarpus arenarius and Salsola subcrassa) 
were all significantly higher under increasing precipitation than under the control in 2011; by 
contrast, no such difference was observed among the precipitation treatments in 2012 (Table 2). 
Based on the community composition, we divided the plant-growing season into an ephemeral 
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dominant phase and a summer-annual dominant phase. Regardless of the precipitation treatment, 
the proportion of plant biomass contributed by ephemerals and summer-annuals did not change, 
as the ephemerals consistently contributed more than 70% of the total community biomass during 
the ephemeral dominant phase, whereas the summer-annuals only contributed to community 
biomass after the ephemerals had died (Fig. 5). Both the P15 and P30 treatments significantly 
increased herbaceous plant biomass in the interplant matrix and beneath shrubs during the 
ephemeral dominant phase (Fig. 5), whereas only the P30 treatment significantly increased 
herbaceous biomass beneath shrubs after the death of the ephemerals (Fig. 5). Increasing 
precipitation significantly increased vegetation coverage by increasing ephemeral coverage 
during the ephemeral dominant phase in 2011 (Fig. 6), yet failed to exert a notable impact in 2012. 
Increasing precipitation apparently did not affect vegetation coverage during the summer-annual 
dominant phase in both years (Fig. 6). 


Table 2 List of plant species within the subplots (1 mx1 m) in the interplant matrix microsite across two years 
(2011 and 2012) in the three precipitation treatments 


Herbaceous density in 2011 Herbaceous density in 2012 

Species/Richness LH (number/m*) (number/m*) 
C P15 P30 C P15 P30 

Erodium oxyrrhynchum E 53x8.0* 67+8.0° 64+3.0° 1241.1" 13x1.0* 18+1.3° 
Alyssum linifolium E 743" 1343° 12+4° 15+3* 2043% 26+5° 
Malcolmia africana E 6x1.0* 6x1.0* 13x5.0* 1+0.6° 240.5° 31.0° 
Leptaleum filifolium E 18.0x12* 12.0x3* 19.0x6* 0.2+0.2° 0.5+0.2* 0.4+0.2* 
Schismus arabicus E 30417" 45423? 42+20° 844" 27+10° 48+17° 
Tragopogon ruber P 0.12x0* 0.40+0.2* 0.16x0.1* 0.20+0.2* 1.00+0* 0.60+0.3* 
Nonea caspica E 2.70x1.6* 0.12x0* 0° 5.00x0* 3.00x1* 3.50x1* 
Alyssum desertorum E 1.30x1* 0° 0.04x0* 1.30x1.0* 0° 0.70+0.2* 
Amberboa turanica E 1.90x0.6* 3.30+0.8" 0.04+0* 1.80x0.8* 4.50+1.5° 4.20+1.0° 
Arnebia guttata E 0? 0° 0.04+0° 0° 0° 0° 
Nepeta micrantha E 0.30+0. 1" 0" 0.04x0* 0? 0.80+0.2* 0° 
Ceratocarpus arenarius A 8+7.0° 1245.9% 1244.0” 16+1.0° 213.0” 55+8.0° 
Salsola subcrassa A 146+95* 177438" 173448” 1.8+0.4* 0° 0.6+0.2* 
Tetracme recurvata E 0° 0° 0° 0.4+0.2* 0.9+0.2* 0.8+0.2* 
Haloxylon ammodendron P 2.6x1.9* 4.0x1.3* 1.6+0.6* 0.8+0.5* 0.8+0.2" 0.8+0.2* 
Limonium bicolor P 0.08x0* 0.16x0* 0.04x0* 0.80+0.6" 0.40x0* 0.50x0* 
Salsola passerina S 8.426.0* 5.9+4.9° 4.2+2.5° 12.0+2.0° 4.0+1.0° 1.0+0.2° 
Descurainia sophia E 0.12x0* 0° 0.04x0* 0° 0° 0° 
Chenopodium glaucum A 0.04x0* 0° 0° 0° 0° 0° 
Euphorbia turczaninowii A 0° 0° 0° 0.40x0.1* 0° 0.16+0° 
Salsola ruthenica A 0.12+0° 0.04+0° 0.04x0* 0? 0? 0.50x0* 
Carex physodes P 0° 1-24" 0° 0° 0° 0.8+0° 
Lactuca undulate A 0.12x0* 0.08x0* 0.08x0* 0° 0° 0° 
Hyalea pulchella A 0.40+0.04* 0.48+0.20* 0.2040. 10# 0.40x0* 0.20x0* 1.50x0.50* 
Seriphidiam santolinum P 0.40x0.2^ 0.12+0* 0.20x0.1* 0° 0.50+0° 0.20+0° 
Koelpina linearis P 0.04+0° 0° 0° 0° 0.80+0° 0° 
Richness 12.0+0.3* 13.4x1.0* 14.0x1.0* 8.8+0.6° 10.0+0.6a 11.8+0.7* 


Note: Different lowercase letters indicate significant difference between treatments within rows for each year at P«0.05 level. LH, A, P, 
and E respectively indicate life history, annual, perennial, and ephemeral plants. 
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Fig. 5 Herbaceous community biomass in response to increasing precipitation in the interplant matrix (a, b) and 
beneath the shrubs (c, d) of Haloxylon ammodendron in 2011 and 2012. The plant community's growing season 
could be divided into an ephemeral dominant phase (grey) and a summer-annual dominant phase (white). Insets in 
each panel show the proportion of biomass of ephemerals (E) and summer-annuals (A) in the two phases under 
the three precipitation treatments, measurements were conducted at the growth peak. C, control; P15, extra 1596 
precipitation addition; P30, extra 30% precipitation addition. 
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Fig. 6 Herbaceous vegetation coverage in response to increasing precipitation at the growth peak of the 
ephemerals and the summer-annuals in 2011 (a) and 2012 (b). C, control; P15, extra 15% precipitation addition; 
P30, extra 30% precipitation addition. 


3.5 Factors affecting the ecosystem WUE 


WUE was positively correlated with GEE and ET (Figs. 7a and b) in 2011, whereas WUE was 
negatively correlated with ET in 2012 (Fig. 7b). No significant relationship was found between 
WUE and plant biomass (Fig. 7c). However, for the interplant matrix, WUE was positively 
correlated with vegetation coverage in 2011 but this relationship was negative and at near 
significance in 2012 (Fig. 7d). 


4 Discussion 


4.1 Ecosystem WUE in response to increasing precipitation 


The WUEuwzz varied from —2.71 to —1.00 g C/kg H20 in our study. Using the eddy covariance 
technique, we measured WUEnre and the measured values ranged from —2.30 to 0.90 g C/kg H20 
in a 2-year study with annual precipitations of 147 and 197 mm in a desert shrub ecosystem 
(Hastings et al., 2005). Using the eddy covariance technique at our study site, Liu et al. (2012) 
measured WUENzz and reported that the values were 0.03 and 0.15 g C/kg H5O under annual 
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precipitations of 127 and 231 mm, respectively. In our study, the negative WUEwzz was resulted 
from the negative values of NEE, suggesting that precipitation-induced ecosystem carbon release 
exceeds that of carbon uptake via herbaceous plant growth. The positive values for NEE in this 
desert measured from the eddy covariance method mainly occurred during the night, largely due 
to the downward carbon absorption by the alkaline soil (Ma et al., 2012). As our carbon exchange 
measurements were taken only in the daytime, the NEE can represent the carbon release occurring 
most of the time in the growing season. 
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Fig. 7 Linear relationships between ecosystem water-use efficiency (WUE) and ecosystem carbon exchange 
(GEE, a), ecosystem evapotranspiration (ET, b), plant biomass (c), and vegetation coverage (d) in the interplant 
matrix and beneath the shrubs of Haloxylon ammodendron (“Beneath H. a”) in 2011 and 2012 


WUEnez responded differently to increasing precipitation in the two study years. In the normal 
precipitation year (2011), increasing precipitation significantly increased WUEnwzz; however, 
increasing precipitation inhibited WUEnwzz in the dry year (2012). The results illustrate that the 
effects of increasing precipitation on ecosystem WUE depend on its stimulation of plant growth 
(Cohen, 1970; Golluscio and Oesterheld, 2007; Liu et al., 2012). In the normal precipitation year 
(2011), 67.9 mm of precipitation fell during the ephemeral growing season, representing 54.3% of 
the precipitation for the entire growing season. Correspondingly, increasing precipitation in 2011 
met the water requirement of the spring ephemerals, leading to stark increases in NEE and 
decreases in ET. In the dry year (2012), when only 20.9 mm of precipitation fell during the 
ephemeral growing season, a large amount of precipitation occurred in the summer-annual 
growing phase; thus the increasing precipitation applied was able to meet the water requirement 
of the summer-annuals. The relative importance of the two functional herb groups can be seen in 
terms of their biomass and coverage: ephemerals contributed to more than 7046 of the herbaceous 
biomass. Therefore, increasing precipitation exerted a stronger stimulation on NEE than on ET, 
leading to the increased WUE seen in 2011. However, increasing precipitation resulted in a 
decreased NEE and similar ET, leading to a decreased WUE in 2012. 

Our study implies that changes in community composition can influence the response of WUE 
to increasing precipitation. Although WUE showed a positive correlation with precipitation 
amounts at the ecosystem level (Varnamkhasti et al., 1995; Ruppert et al., 2012), increasing 
precipitation may exert a different effect on different vegetation types (Huxman et al., 2004). In 
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addition, the seasonal pattern of the precipitation regime should be considered when predicting 
precipitation effects on ecosystem WUE. 


4.2 Leaf WUE and ecosystem WUE respond differently to increasing precipitation 


Leaf WUE (the ratio of leaf CO» assimilation rate to the transpiration rate) has been widely 
studied through direct measurement of instantaneous gas exchange rates at the leaf level, often 
with portable equipment and a carbon isotopic measurement used for longer time periods (Toft et 
al., 1989; Wang et al., 2005; Golluscio and Oesterheld, 2007; Wu et al., 2016). Leaf WUE is 
controlled principally by leaf stomatal dynamics. Although leaf WUE varies at daily and seasonal 
scales, the values measured at certain times (i.e., at mid-morning) are used to represent the 
whole-day WUE. To date, we have not reached consensus yet as for how leaf WUE responds to 
increasing water availability in arid and semi-arid regions (Cohen, 1970; Ehleringer and Cooper, 
1988; Toft et al., 1989; Lin et al., 1996; Golluscio and Oesterheld, 2007). In our study, WUEL 
slightly decreased with the increasing precipitation, regardless of plant species, which suggests a 
stronger stimulation of leaf transpiration than photosynthesis. 

To be partial supportive to our first hypothesis, WUEmwzz responded differently from WUE, to 
the increasing precipitation treatments in our study. WUENwzz increased with increasing 
precipitation, while WUE, slightly decreased in 2011. Unfortunately, little is known about how 
leaf-scale variation in WUE influences landscape-scale variation in ecosystem WUE (Linderson 
et al., 2012). On one hand, due to a large spatial variation of leaf WUE within and among 
different vegetation canopies, it is difficult to reliably extrapolate leaf WUE to the whole-plant 
WUE and ecosystem WUE (Niu et al., 2011; Medlyn et al., 2017). On the other hand, ecosystem 
respiration and soil evaporation are two important factors controlling ecosystem WUE, but how 
they responded to increasing precipitation in our study remained unknown. The relationship 
between leaf WUE and ecosystem WUE was reportedly different from some studies to other 
studies and the difference was probably caused by different factors considered in different studies. 
In brief, the previous studies show that leaf WUE hardly explained the variation of ecosystem 
WUE (Niu et al., 2011; Medlyn et al., 2017). The response of ecosystem WUE to increasing 
precipitation is strongly dependent on the effects of increased water availability on the recipient 
herbaceous community, and water availability determines ecosystem NEE and ET in our study. 
The increased NEE in 2011 was clearly related to the strong increase in plant biomass. 
Community biomass under the extra 30%-increasing precipitation (P30) was nearly 2.5 times that 
of the control, and this increased plant biomass contributed to the greater increase in GEE than in 
ER. Meanwhile, the increased vegetation cover under increasing precipitation also decreased the 
soil evaporation. Thus, the increased biomass and decreased soil evaporation together generated 
the positive response of ecosystem WUE to the increasing precipitation. In contrast, an increased 
GEE was insufficient to offset ER and to decrease soil evaporation under the increasing 
precipitation in 2012; thus, ecosystem WUE was not altered by increasing precipitation. Our 
results indicate the ecosystem WUE response to increasing precipitation is related not only to the 
precipitation amount but also to the precipitation seasonality and plant functional groups. 


4.3 Responses of ecosystem WUE to increasing precipitation under two contrasting 
microsites 


There was no apparent difference in WUENwzz between interplant matrix and beneath the shrubs of 
H. ammodendron in the normal precipitation (2011) and dry (2012) years. This result does not 
support our second hypothesis: ecosystem WUE is higher in the interplant soil microsite than that 
beneath shrubs in the normal precipitation year, and ecosystem WUE is similar in the dry year 
both in the interplant matrix and the beneath shrubs. No differences in ecosystem WUEzzx 
between two habitats are attributable to the similar plant biomass and vegetation cover. Previous 
studies have shown that an alteration of the dominant species can cause marked changes in 
ecosystem WUE, which is not only related to species-specific photosynthetic traits but also to the 
alteration in soil respiration and soil H2O fluxes induced by a changed canopy structure 
(Golluscio and Oesterheld, 2007; Niu et al., 2011; Eamus et al., 2013). Community composition 
in the interplant matrix microsite was different from that found beneath shrubs and was largely 
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driven by rare and accompanying species, whereas the dominant species were the same or had 
similar functional traits both in the interplant matrix and the beneath shrubs. For example, in the 
ephemeral growing phase, E. oxyrrhynchum comprised most of the community biomass of the 
interplant microsite, while S. arabicus and C. glaucum were the major contributors to the 
vegetation beneath shrubs. However, in the summer-annual growing phase, Salsola subcrassa and 
S. passerine were the major contributors, respectively, to the community biomass of the interplant 
matrix and beneath shrubs. NEE and ET showed no discernible difference between the two 
microsites in the normal precipitation year (2011), leading to similar ecosystem WUEs. In the dry 
year (2012), NEE and ET showed synchronous increases for the interplant matrix, giving rise to 
similar ecosystem WUEs between the two microsites. The dominant species were always similar 
between the two habitats, thus community structure and coverage were also similar. By comparing 
the results of ecosystem WUE between typical microsites in a desert ecosystem, this study 
revealed that WUE was similar between the interplant matrix and beneath shrubs in this desert in 
China. 


5 Conclusions 


In this study, we found that leaf WUE and ecosystem WUE showed inconsistent responses to 
increasing precipitation. Leaf WUE consistently decreased as precipitation increased. By 
contrast, the responses of the ecosystem WUE to increasing precipitation are different in 
different precipitation regimes: increasing by 33.9% in the wet year (i.e., the normal 
precipitation years) and decreasing by 4.196 in the dry year when the precipitation was about 3096 
less than that in the wet year. Increasing precipitation increased plant biomass, but it did not shift 
the relative significance of ephemerals in this vegetation community, thus increasing ecosystem 
WUE in 2011. By contrast, increasing precipitation exerted no significant effects on ecosystem 
WUE in 2012 because plant biomass could not respond to it in a dry year. Our field study shows 
that the effects of increasing precipitation on ecosystem WUE depended heavily on the 
stimulating effects of the dominant species in the plant community. Hence, altering the 
community composition through changing precipitation regime could predictably shift 
ecosystem WUE in this desert ecosystem. 
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